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ABSTRACT 


(  This  report  presents  some  of  the  results  of  the  experimental 
and  theoretical  studies  on  the  fatigue  crack  propagation  in  cy- 

fS  J*  ^  f'  ***  ^  ■■  

ft  T  "  7  ^  ^  ^  \ 

lindrical  shells^subJeTteT  to  fTu“^^^  internal  pressure. J 

Most  of  the  fatigue  tests  are  performed  on  6063-T6  aluminum 
alloys.  Some  2024-T3  aluminum  and  mild  steel  specimens  are  also 
tested.  Plexiglas  specimens  are  used  for  static  experiments. 


[Tn  the  analysis  of  the  experimental  results,  the  stress  in¬ 
tensity  factor  is  used  as  the  correlation  parameter.  A  modified 
crack  propagation  model  is  developed  to  take  into  account  the 
effect  of  bending  stresses,  which  are  superimposed  on  membrane 
loads  around  the  crack  in  shell s .  |  To  evaluate  the  stress  inten¬ 
sity  factor,  the  integral  equations  obtained  by  Folias  are  solved 


numerically  after  separating  the  singularities.  Partly  to  ver¬ 
ify  the  theoretical  results,  strain  measurements  were  made 
around  the  crack  tip  on  the  outside  of  the  aluminum  and  inside 
and  outside  of  the  plexiglas  cylinders.  The  results  indicate' 
that  over  the  range  of  crack  propagation  data  covering  10  ^  to 
10"^  in/cycle,  the  model  adopted  in  this  study  seems  to  be 
highly  satisfactory. 
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INTRODUCTION 


The  problem  of  fatigue  crack  propagation  in  thin  plate  and 
shell  structures  is  of  considerable  technological  importance, 
as  it  relates  directly  to  the  fail-safe  design  of  aerospace  ve¬ 
hicles  and  some  ship  structures.  In  such  structures,  because 
of  the  unavoidable  existence  of  stress  rai sers  ,  (e . g .  ,  holes, 
various  forms  of  joining,  material  imperfections),  generally  the 
fatigue  crack  nucleates  after  relatively  small  number  of  load 
cycles  compared  to  the  total  useful  life  of  the  structure. 

Hence,  fatigue  crack  propagation  studies  become  necessary  for 
rational  design  as  well  as  maintenance  procedures. 

The  main  objective  of  the  present  study  is  to  investigate 
the  effect  of  shell  curvature,  and  to  a  lesser  extent,  that 
of  bending  stresses  super imposed  on  the  membrane  loads  on  the 
crack  propagation  rate  in  thin-walled  structures.  In  published 
literature  there  is  a  great  wealth  of  data  as  well  as  the  results 
of  studies  on  the  crack  propagation  in  a  large  variety  of  struc¬ 
tural  materials  subjected  to  repeated  uniaxial  tension.  Thus, 
it  becomes  desirable  to  develop  a  method  by  which  the  fatigue 
crack  growth  characteristics  of  thin-walled  structures  may  be 
predicted  from  the  fatigue  information  on  flat  plates  of  the 
same  material  subjected  to  uniaxial  extension. 


Experience  with  the  fatigue  crack  propagation  in  flat 
plates  indicate  that  from  an  engineering  view  point  the  stress 
intensity  factor  for  the  propagating  crack  is  a  very  simple 
and  effective  parameter  in  analyzing  and  correlating  the 
results  of  the  fatigue  tests.  This  will  also  be  the  view 
point  which  will  be  adopted  in  this  study,  particularly  in 
the  modification  of  the  crack  propagation  model  to  take  into 
account  the  effect  of  the  combined  loads.  It  appears  then 
that  there  are  two  main  problems  which  should  be  considered. 
One  is  the  evaluation  of  the  stress  intensity  factor  through 
the  elastostatic  analysis  of  an  internally  pressurized  cylin¬ 
drical  shell  containing  a  longitudinal  crack,  and  the  other 
is  the  experimental  verification  of  the  crack  growth  model 
by  using  the  results  of  the  elastic  analysis  and  the  fatigue 
tests . 

For  the  stress  intensity  factor  an  assymptotic  solution 
is  given  by  Folias  [1,2,3].  However  this  solution  is 
reliable  only  for  very  small  crack  lengths  and  a  more  refined 
solution  is  needed  for  the  crack  dimensions  or,  more  precisely 
for  the  shell  parameters  which  may  be  encountered  in  practice. 
To  obtain  such  a  solution  the  integral  equations  governing 
the  problem  is  solved  numerically  after  separating  the 
si ngul ari ties. 

Most  of  the  fatigue  crack  propagation  data  in  this  study 
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are  obtained  by  using  6063-T6  aluminum  tubes.  In  the 
subsequent  series  of  tests  6061-T3  tubes  and  plates  of  the 
same  thickness  are  being  used. 
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FATIGUE  CRACK  PROPAGATION  MODEL 


Almost  all  of  the  existing  quantitative  fatigue  crack 
propagation  models  deal  with  the  thin  plates  with  straight, 
through  cracks,  which  are  subjected  to  uniaxial  repeated 
extensional  loads  [4-16],  A  summary  of  various  models  as 
well  as  the  description  of  a  crack  propagation  model  based 
on  plastic  deformations  around  the  propagating  crack  were 
given  in  an  earlier  report  [17].  In  this  model  the  crack 
propagation  rate  is  expressed  as 


da  ^  »  ^“1 
HTT  ^  Pmax 


(1) 


where  a  is  the  half  crack  length,  n  is  the  number  of  cycles 
A,a.|,  a2  are  positive  constants,  ^r  respectively, 

the  maximum  and  range  values  of  the  plastic  zone  size  around 
the  crack  tip. 

In  flat  plates,  a  reliable  estimate  of  the  plastic  zone 
size  may  be  obtained  by  using  Dugdale's  technique  [18], 
which  gives 

p  =  a[sec(|| — )  -  1]  (2) 

ys 

where  a“  is  the  uniaxial  stress  perpendicular  to  the  plane  of 

and  away  from  the  crack  and  a  is  the  yield  stress. 

y  s 
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Similarly,  in  the  case  of  cylindrical  bending  the  fatigue 
crack  propagation  rate  may  be  expressed  as  [15,16,17] 

=  B(l+fi)^“l 

where  o.  and  are  given  by  (4)  with  a“  as  the  bending  stress 
on  the  surface  of  the  plate  and  y  has  a  theoretical  value  of 

0.5. 

Another  model,  which  accounts  for  the  effect  of  mean 
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stress  as  well  as  stress  range,  is  that  given  by  Forman  and 
others  [19] 

00 

da/dn  =  C(K^)'"/[(1  -  -  2K^]  (7) 

°mi  n 

where  C  and  m  are  material  constants  and  is  the  critical 
stress  intensity  factor.  The  crack  growth  models  used  in 
this  report  to  express  the  crack  growth  rate  in  cylindrical 
shells  will  be  of  the  same  form  as  (5)  and  (6),  subject  to 
the  limitation  that  small  scale  yielding  exists  at  the  crack 
tip. 

In  the  neighborhood  of  the  fatigue  crack  if  the  plate  is 
subjected  to  a  combination  extension  and  bending,  because  of 
linearity,  the  resulting  stress  field  would  be  the  direct 
addition  of  the  two  separate  fields.  This  means  that  the 
range  component  of  the  resultant  stress  intensity  factor 
would  be  +  K^j^.  It  is  obvious  that  using  this  value 

in  the  crack  propagation  models  expressed  in  terms  of  K, 
e.g.,  (5)  or  (7),  would  give  erroneous  results,  since  it 
would  not  reduce  to  (6)  as  ,  the  extensional  component  of 
the  stress  intensity  factor  goes  to  zero.  Thus,  one  way  of 
superimposing  the  effects  of  extension  and  bending  in  fatigue 
crack  propagation  is  suggested  by  (5)  and  (6),  namely, 
replacing  in  crack  growth  models  by  giving 
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^  =  B(I  +  0)2“1  +  YK^b)^<''l*“2> 


(8) 


In  (8)  it  is  assumed  that  mean-to-range  stress  ratios,  i.e., 
the  values  of  fi  in  bending  and  extension  are  the  same.  This 
is  the  case  in  shells.  If  values  in  bending  and  extension 
are  not  equal,  (8)  would  have  to  be  modified  accordingly. 

Note  that  (8)  may  be  obtained  by  assuming  that  for  small 
scale  yielding,  the  plastic  zone  size  in  the  neighborhood  of 
a  crack  in  a  plate  subjected  to  combined  plane  extension  and 
cylindrical  bending  loads  may  be  estimated  as 

P  -  1  (9) 

Due  to  the  nonlinear  nature  of  the  phenomenon,  methods  of 
superposition  of  the  effects  of  combined  loads  based  on  the 
addition  of  plastic  zone  sizes  or  crack  growth  rates  are  not 
justifiable. 


ELASTIC  SOLUTION  FOR  CRACKED 
CYLINDRICAL  SHELLS 

The  details  of  the  elastic  solution  for  a  cylindrical 
shell  containing  a  longitudinal  crack  are  given  in  Appendix  I. 
This  is  an  extension  of  the  solution  given  by  Folias  [2].  It 
is  restricted  to  shallow  shell  theory  of  isotropic,  homo¬ 
geneous,  elastic  circular  cylinders  of  constant  thickness. 
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In  this  case,  the  system  of  differential  equations  derived  by 
Marguerre  [20]  takes  the  form 


Eh 

R  JyJ 


= 


0 


(10) 


v'+W 


i_  a^F  _  q(X,Y) 
RD  JyJ  '  D 


where  W  is  the  displacement  along  the  Z  axis,  F  is  the  stress 
function,  X,  Y,  Z  are  the  rectangular  coordinates  shown  in 
Figure  1,  E  is  the  modulus  of  elasticity,  h  is  the  shell  wall 
thickness,  R  is  the  shell  radius,  q  is  the  internal  pressure, 
and  D  is  the  flexural  rigidity  defined  by 


D  =  Eh3/[12(l-v2)] 


The  bending  moment  components  M  ,  M  ,  and  M  are  given 
by 


Mx  = 


D[ 


92W 

3X^ 


32Wn 

V 


M  = 

y 


nrB^W 

DL^ 


32W. 

ixf-' 


(11) 


M 


xy 


=  -  D(l-v) 


92W 

3X3Y 
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where  v  is  Poisson's  ratio.  The  membrane  forces  N^,  N^,  and 
are  given  as 


N 


X 


N 


y 


a^F 

J77 


(12) 


a^F 

xy  "  aXaY 


and  the  transverse  shear  forces  and  are  given  as 


Qx  =  -  D  (v^W) 


(13) 


Q 


y 


D  Iy  (v2W) 


In  view  of  equations  (11),  the  bending  stress  components 
become 


®xb 
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(14) 


^xyb 


2GZ 


a^w 

axaY 


9 


and  from  equations  (12)  the  extensional  stress  components  are 
obtained  as 


1  92f 
^xe  “  h 

_  1  a^F 
°ye  h  aX^ 


(15) 


1  a^F 

^xye  "  ■  h  aXaY 

Introducing  the  dimensionless  coordinates 


(16) 


where  a  is  the  half  crack  length,  the  homogeneous  parts  of  (10) 
become 


Eha2  a2W  , 
R  ax2 


v‘+F  =  0 


(17) 


ai  a2F  _ 
RD  ax2 


0 


/The  solution  for  the  internally  pressurized  cracked  cylinder 
is  obtained  by  first  evaluating  the  stresses  in  a  cylinder  with¬ 
out  a  crack,  obtaining  the  bending  moments,  membrane  forces  and 
transverse  shear  1 oads at  the  location  of  the  crack,  and  then  ap¬ 
plying  equal  and  opposite  of  these  moments  and  loads  on  the  sur-  ) 
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face  of  the  crack  in  a  cylinder  with  no  other  external  loads. 
Superposition  of  these  two  results  gives  the  desired  sol uti on . J 


Since  the  first  problem  has  no  contribution  to  the  singularity 
at  the  crack  tip,  and  since  we  are  interested  in  evaluating  the 
stress  intensity  factor  only,  the  main  problem  is  the  solution 


of  the  homogeneous  system  (17)  with  nonvanishing  loads  at  the 
crack  surfaces. 


In  solving  (17),  additional  requirements  to  be  met  are  that 
the  displacement  function  W  and  stress  function  F  with  their 
first  derivatives  must  be  finite  far  away  from  the  crack  and  the 
stresses  and  displacement  for  y  =  0  and  |x|>l  must  be  continuous 
i.e.,  F  and  W  and  all  their  partial  derivatives  must  be  contin¬ 
uous  for  all  X  and  y  except  for  points  on  the  crack  surface, 


1  i  m 
|y  1^0 


ay  ay 


0 


,  (n  =  0,1, 2, 3)  (18) 


lim  [^(f'^)  -  )]  =  0 

I  y  I  ^0  ay  ay 


For  the  symmetric  case,  N  and  the  effective  transverse 

3M  a-t  y=-o 

shear  resultant,  a>"e  zero^and  the  functions  W 

and  F  may  be  represented  in  terms  of  Fourier  cosine  integrals. 
Then  the  boundary  conditions  and  the  continuity  requirements 
lead  to  a  system  of  dual  integral  equations.  These  dual  inte- 
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gral  equations  can  be  reduced  to  singular  integral  equations 
with  kernels  L-i,  and  L^.  These  singular  integral  equa¬ 

tions  have  the  form 


1  2vH  a2  1 

/  J_  dx  =  /  [u,{5)  L,(5,x)  +  Up(c)  Lp(?,x)]  d? 

-1  i/EhD  -1  '  '  ^ 


(19) 


/  a2TTMy  §  L3(?,x)  +  02(5)  L^Ce.x)]  d? 

n  Dm 

where  =  -^  and  (n^  and  m^  are  chosen  as  constants). 

The  solutions  u.j(c)  and  U2(c)  may  then  be  expressed  as 


U■^{K)  =  [A^  +  x2A^(l  ^2)  +  ..  J.  |5|<1 

(20) 

1^2(5)  "  [Bq  +  x2B.|(1-c2)  +  ...];  |5|<1 

where  Ap  and  Bp  are  complex  constants  and 
X  =  [12(l-v2)]^/4a/(Rh)^/^ 
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A  and  Bp  are  generally  functions  of  x.  Substituting  (20)  into 
(19)  and  writing  the  kernels  in  series  form,  good  for  the  range 
0<x<2,  the  resulting  integrals  can  be  solved  in  closed  form. 

The  resulting  values  of  the  integrals  are  polynomials  in  odd 
powers  of  x  of  the  following  form 


V  r  ^2k+l 

“  I 

i/EhD  0 

(21) 

-  2 Tim  X  =  I  D|^x^ 

0 

where  Cj^  and  D|^  are  functions  of  X,  Ap  and  Bp.  By  comparing 
coefficients  of  x^ ,  a  system  of  simultaneous  equations  can  be 
obtained  and  solved  for  Ap  and  Bp. 

The  strength  of  the  singularities  at  the  crack  tips  depend 
on  the  constants  A^  and  B^  only,  and  the  stresses  in  a  small 
neighborhood  of  the  crack  tip  may  be  expressed  as 
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(22) 
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xye 
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(23) 


where  P||^  and  P^  are  functions  of  and  is  the  stress  in¬ 

tensity  factor  for  a  plate  under  extensional  loads  (K^  =  a“/a) , 
and  r  is  the  distance  from  the  crack  tip.  The  values  +P|!^  and 
-P|Ij  are  used  for  stresses  on  the  inner  and  outer  surfaces  of  the 
shell,  respectively.  As  a->0  (i.e.,  R^«>) ,  P^->0  and  P^-^l  and  the 
stresses  of  a  flat  plate  in  extension  and  bending  are  recovered. 
Thus  the  stresses  in  the  crack  tip  region  of  the  shell  are  ex¬ 
pressed  in  terms  of  stresses  in  a  flat  sheet.  The  stress  inten¬ 
sity  factors  for  a  shell  with  a  meridional  crack  due  to  bending 
and  extension  become,  respectively 


Kbs  =  [-P^//(l-v'^)/3r 


(24) 


K 


es 


e  e 


(25) 
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Note  that  because  of  the  Kirchhoff  boundary  conditions, 
the  bending  shear  stress  does  not  vanish  at  the  free  edge  of 
the  crack. 

Since  is  the  stress  intensity  factor  of  a  plate  under 
extensional  loading,  the  combined  stress  intensity  factor,  , 
for  a  cylinder  containing  a  meridional  crack  can  be  written  as 


for 

the 

outer 

surface 

(26) 

for 

the 

inner 

surface 

(27) 

where  "  P^//(l-v^)/3 

Therefore,  the  values  Pg  and  combine  to  form  the  stress 
intensity  ratio  of  a  shell  under  internal  pressure  to  a  plate 
under  extension,  for  a  particular  material. 

Computed  values  of  and  Pg  are  given  in  Table  I  for 
values  of  A  from  0.2  to  2.2  in  increments  of  0.2  and  for  num¬ 
ber  of  terms  of  A  and  B  in  the  series  expansion  (20).  As 

r  • 

could  be  expected,  the  convergence  becomes  slower  as  x  is  in¬ 
creased.  For  x>2,  the  convergence  is  not  expected  to  be  good 
because  of  the  limits  imposed  on  x  in  this  particular  series 
form  of  the  kernels.  A  minimum  accuracy  of  three  digits  was 
achieved  for  all  Pg  for  .2_<X£2.2.  This  was  not  the  case  for 
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was  con- 


P|I|.  But,  since  Pg>>|P|l,|.  the  accuracy  achieved  for  Pjl^ 

sidered  sufficient.  Figures  3  and  4  are  plots  of  P  and  P ' : 

e  b 

and  P|^  for  v  =  1/3  obtained  from  the  asymptotic  solution  [2], 
a  two  term  and  an  eight  term  expansion  in  (20).  In  the  asymp¬ 
totic  solution,  orders  of  a'+  and  greater  were  neglected  in  the 
kernels  and  the  solution  for  and  B^.  Four  simultaneous 
equations  were  used  in  solving  for  and  B^.  Figures  3  and 
4  compare  this  result  with  the  two  and  eight  term  solutions, 
which  consisted  of  the  solution  of  four  and  sixteen  simulta¬ 
neous  equations.  In  the  latter,  the  highest  degree  term  in  a 
which  is  neglected  is  a^s.  jhe  eight  term  solution  shows  that 
the  sign  change  in  the  bending  stress  intensity  ratio  P^  (Fig¬ 
ure  4)  does  not  take  place  in  the  interval  under  consideration. 
This  means  that  for  this  a  range,  the  maximum  stress  in  the 
crack  region  is  on  the  outside  of  the  shell. 

A  least  squares  curve  fitting  program  was  used  to  express 
P|I^  and  Pg  as  polynomials  in  a  which  are  given  below. 

P||j  =  -  5.6207  X  10"^  +  .0016268a  -  .14611a2 

+  .25010A3  -  .2172^**  +  .09771  1a5  (28) 

-  .017632a6 
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Pg  =  1.0004  -  .062432a  +  .67419x2  -  1.0367x3 

+  1.0184x‘*  -  .49949x3  +  .094455x6  (29) 


A  comment  is  in  order  concerning  the  use  of  the  Kirchhoff 
theory  in  the  above  solution,  according  to  which  only  the  ef¬ 
fective  transverse  shear  rather  than  Q  and  M  is  required  to 

y 

vanish  along  the  crack.  While  the  stress  distribution  outside 
[  the  local  region  should  be  accurate,  it  might  be  expected  that 

near  the  crack  boundary  the  same  type  of  difference  exists  as 
that  found  by  Feridun  [21],  Knowles  and  Wang  [22]  and  Hartranft 
[23]  in  comparing  the  stress  distributions  due  to  bending  using 
Kirchhoff's  and  Reissner's  theories.  In  this  case,  the  order 
of  the  stress  singularity  was  the  same  for  both  theories  but 
the  stress  distribution  in  the  crack  region  was  different, 
Reissner's  theory,  which  accounts  for  all  the  physical  boundary 
conditions,  giving  the  same  angular  dependence  as  the  plane 
stress  solution. 

The  important  fact  is  that  the  ratio  of  the  stress  inten¬ 
sity  of  a  cylinder  to  that  of  a  plate,  +  Pj^,  is  the  desired 
result  for  use  in  the  crack  growth  model  and  that  the  ratio 
will  be  the  same  no  matter  which  one  of  the  two  theories  is 
>se^  to  establish  the  ratio  in  bending.  In  addition,  the  ef- 


17 


feet  of  bending  is  small  for  0£X<2.2  where  P,^  varies  from  0 
to  10%  that  of  Pg. 
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EXPERIMENTAL  PROCEDURE 


The  experiments  were  conducted  to  determine  the  crack 
propagation  rates  and  the  stresses  caused  by  the  crack  in 
thin-walled,  cylindrical  shells  under  different  loading  con¬ 
ditions.  The  loading  employed  in  all  tests  produced  various 
applied  mean  and  range  stresses  from  one-eighth  to  one-quarter 
of  the  yield  stress. 


Description  of  Test  Specimens.  -  Three  materials  were 
tested  with  varying  degrees  of  success  -  2024-T3  aluminum, 
6063-T6  aluminum,  and  mild  steel.  The  best  results  were  ob¬ 
tained  from  the  tests  of  the  6063-T6  aluminum.  Thin-walled 


2024-T3  bare  aluminum  tubing  could  not  be  purchased  commercial¬ 
ly.  So,  these  specimens  were  rolled  into  a  cylinder  from  a 
flat  sheet.  The  rolled  cylinder  was  8  5/8  inches  in  diameter, 
27  inches  in  length  and  0.050  of  an  inch  in  wall  thickness. 


The  mating  edges  as  well  as  the  ends  of  the  shell  were  clamped 

to  a  mount  to  form  the  pressure  chamber.  pThe/only  thin-walled 

tubing  that  was  commerically  available  was/6063-T6  extruded 

i  A 

aluminum  cyl i nder^.  J  Eight  test  specimens  were  cut  from  one 
20  foot  long  tube.  The  mechanical  properties  and  dimensions 
of  the  6063-T6  aluminum  specimens  are  given  in  Tables  II  and 


III.  The  steel  shells  were  rolled  from  a  0.050  inch  thick  flat 
plate  into  an  8  inch  diameter  cylinder  and  then  seam  welded, 
fjence,  only  the  ends  of  the  6063-T6  aluminum  and  the  steel 
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shells  needed  to  be  clamped  to  a  mount  to  form  the  pressure 
chamber. 

Experimental  equipment.  -  The  experimental  set-up  and  a 
schematic  diagram  of  the  hydraulic  system  are  shown  in  Figures 
5  and  6,  respectively. 

The  pressurization  control  system  basically  consisted  of 
the  control  panel,  control  manifold,  and  hydraulic  power  sup¬ 
ply.  It  maintained  any  desired  pressure  loading  characteristics 
in  the  specimen  by  means  of  a  simple  proportional  feed  back  con¬ 
trol  system.  The  function  generator,  located  in  the  control 
panel,  determined  the  shape  of  the  pressure  variation  curve  and 
the  rate  and  character  of  the  pressure  loading.  The  mean  and 
range  pressures  were  established  by  the  control  panel  s  set 
point  and  span  controls. 

A  50  power  microscope,  mounted  on  a  4  inch  micrometer  stage, 
was  used  for  crack  readings.  It  had  a  smallest  division  reading 
of  5x10'^  inches,  which  was  actually  more  accurate  than  needed 
since  the  crack  tip  was  not  definable  to  more  than  10  inches. 

The  strain  indicator,  used  for  static  strain  measurements, 
was  accurate  to  0.1%  of  the  reading  or  5  microinches/inch, 
whichever  was  greater.  The  oscilloscope,  used  for  dynamic 
strain  measurements  had  a  smallest  division  reading  of  4  micro¬ 
inches/inch. 
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Speci men^^preparati on .  -  A  crack  was  initiated  in  the  shell 
by  first  drilling  a  0.020  inch  diameter  hole  at  the  point 
representing  the  center  of  the  crack,  and  then  making  a  cut 
that  was  symmetrical  about  the  hole  with  a  0.011  inch  wide 
jewelers,  saw.  A  notch  was  then  formed  at  each  end  of  the  saw 
cut  by  "sawing"  with  a  razor  blade.  The  minimum  total  length 
of  the  manufactured  crack  was  0.200  inches.  I  Next,  the  speci- 
men  was  sanded  and  polished  in  the  crack  region  to  remove  the 
surface  flaws  and  make  the  crack  tip  visible. 


Since  the  shells  were  pressurized  with  oil,  the  crack  had 
to  be  sealed  with  a  seal  that  had  little  or  no  effect  on  the 
stress  distribution  of  the  shell  in  the  crack  region.  This 
was  accomplished  with  a  steel  shim  and  vinyl  patch.  1  The  shim 


was  .003  inch  thick  and  about  2  inches  by  3  inches  in  size,  was 
placed  behind  the  crack  on  the  inner  portion  of  the  shell,  and 
was  coated  with  a  thin  layer  of  graphite  to  eliminate  any  pos¬ 
sible  effect  of  friction.  The  vinyl  patch  was  .008  inch  thick 
and  about  6  inches  square,  and  was  cemented  to  the  shell  and 
to  the  center  of  the  shim.  Thus  the  shim  was  prevented  from 
shifting  from  behind  the  crack  while  still  free  to  slide  on 
the  shell  and  to  conform  to  the  shell  contour  in  the  crack 
region. 


Specimen  mounting.  -  The  shell  was  mounted  on  a  machined 
steel  pipe.  Figure  7  is  a  cross-sectional  drawing  of  the 
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shell  and  mount  for  the  6063-T6  aluminum  and  the  steel  shells. 
The  mount  for  the  2024-T3  aluminum  was  basically  the  same  ex¬ 
cept  for  additional  holes  for  clamping  the  edges  of  the  shell. 

A  layer  of  cement  was  applied  to  the  oil  pressure  side  of  all 
contact  surfaces  of  the  shell,  mount,  and  spacer  to  keep  oil 
leakage  to  a  workable  minimum.  The  rest  of  the  contact  sur¬ 
faces  were  bare  so  that  when  pressure  was  applied,  the  surfaces 
would  not  slip  and  the  strain  in  the  x  direction  would  be  zero. 

The  end  clamps  were  1/2  inch  wider  than  the  inside  contact 
surface.  This  1/2  inch  overhang  had  a  two  degree  chamber  so 
that,  when  the  shell  expanded  under  pressure,  there  would  be 
no  high  stress  concentration  in  the  shell  at  its  ends.  A  soft 
aluminum  (1100-0)  split  ring  was  used  between  the  outer  surface 
of  the  shell  and  end  clamps  to  uniformly  distribute  the  clamping 
force.  After  the  shell  was  clamped  to  the  mount,  strain  gages 
were  cemented  to  the  shell's  outer  surface.  The  gages  were 
BLH  type  A-5  with  nominal  resistance  of  120  ohms  and  gage  fac¬ 
tor  of  2. 


Experimental  procedure.  -  The  test  results  of  the  2024-T3 
aluminum  did  not  give  reproducible  crack  growth  results,  either 
due  to  the  necessity  of  clamping  down  the  edges  as  well  as  the 
ends  of  the  cylinder,  or  due  to  a  large  geometrical  deviation 
from  a  true  cylinder.  However ,|^eproducti bl e  crack  growth  re¬ 
sults  were  obtained  for  the  6063-T6  al umi num.  |  No  conclusions 
could  be  reached  on  the  steel  because  of  the  small  amount  of 


22 


data  collected.  Therefore,  this  report  will  use  only  the  re¬ 
sults  from  the  6063-T6  aluminum  unless  stated  otherwise. 

The  pressure  was  calibrated  with  a  dead  weight  tester  to 
insure  the  proper  prediction  of  the  internal  pressure  of  the 
shell  from  the  control  settings.  The  output  was  found  to  be 
linear  with  the  control  settings.  Maximum  variation  from  lin¬ 
earity  was  found  to  be  10  microinches/inch.  Since  the  readings 
were  greater  or  less  than  a  linear  relationship,  it  was  felt 
that  most  of  this  was  caused  by  small  deviations  in  manually 
setting  the  controls  and  that  this  error  would  be  self-compen¬ 
sating. 


The  pressurization  unit  was  found  to  give  a  null  or  zero 
shift  in  the  pressure  curve.  Therefore,  constant  monitoring 
of  the  unit  to  maintain  the  zero  position  was  required  during 
the  crack  growth  tests. 

/^The  mean  stress  and  stress  range  varied  from  shell  to 
shell  but  were  held  constant  for  a  given  shell.  The  stress 
varied  sinusoidally  with  a  cyclic  speed  of  8/9  cps  and  a  var¬ 
iation  in  peak  to  peak  stress  of  less  than  1%.  No  tests  were 
conducted  with  a  maximum  hoop  stress,  a”  (applied  stress), 
greater  than  one-half  the  yield  stress.  /  The  hoop  stress  is 
found  from  the  relation  , 


(30) 
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where  q  is  pressure,  R  is  internal  radius,  and  h  is  wall  thick¬ 
ness  . 

For  a  cylinder  with  fixed  ends,  i.e.,  e”  =  0,  we  have  the 

J\ 

following  additional  relations 


00 


(31) 


“  _  (l-v2)  qR 

''y  ■  E  h 


(32) 


where  v  is  Poisson's  ratio. 

Strain  gages  were  added  to  the  outer  surface  of  the  shell 
to  measure  the  hoop  strains.  The  gages  were  placed  along  the 
prolongation  of  the  crack  where  they  would  be  unaffected  by 
the  perturbed  strains  due  to  the  end  clamps  or  the  crack. 
Strain  measurements  were  then  calibrated  with  the  internal 
pressure  through  equation  (32)  and  the  calibration  curve  is 
shown  in  Figure  8.  The  maximum  static  pressure  used  in  cali¬ 
bration  was  not  greater  than  the  maximum  test  pressure  for 
that  shell.  This  prevented  any  unwarranted  strain  hardening 
in  the  crack  tip  area. 

The  measured  strain  gage  readings  gave  differences  from 
the  theoretical  of  as  much  as  11%  for  the  6063-T6  aluminum, 
and  25%  for  the  2024-T3  aluminum.  To  find  an  explanation  for 
this,  an  investigation  of  the  change  in  geometry  of  the  outer 
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surface  of  the  shell  due  to  a  change  in  internal  pressure  was 
conducted. 

The  outer  surface  displacement,  aR,  for  a  perfect  cylinder 
with  uniform  wall  thickness  should  be 

AR  =  R  (33) 

Substituting  from  (32)  into  (33),  we  find 

.p_l-v^R^  /o/|\ 

AR  =  -g —  q  (34) 

The  outer  surface  displacement  of  the  test  cylinders, 
taken  with  a  dial  indicator,  did  not  agree  with  the  results 
predicted  by  (34).  The  readings  were  taken  at  a  position  away 
from  the  influence  of  the  crack  and  end  clamps.  The  readings 
were  very  erratic  as  the  pressure  was  increased.  In  some 
places  the  displacement  would  decrease  while  at  other  places 
it  would  increase  more  than  predicted  by  (34).  This  was  due 
to  the  fact  that  initially  the  cylinder  was  not  perfectly  round 
and  straight. 

Even  though  the  strain  gage  readings  were  different  from 
the  theoretical,  they  were  useful  in  insuring  constant  range 
and  mean  stress.  A  minimum  of  two  gages  at  different  locations 
on  the  shell  surface  were  used  to  measure  the  hoop  strains 
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throughout  the  tests  in  case  one  gage  failed.  None  of  the 
gages  failed  during  the  tests  and  the  readings  were  found  to 
be  very  reproducible. 

Foil  strain  gages  with  1/16  inch  square  grids  were  used 
to  determine  the  strains  ahead  of  the  crack  tip  in  plexiglas 
G  and  in  shells  #12,  #14  and  #18  of  6063-T6  aluminum.  Each 
gage  gave  the  average  strain  of  the  surface  covered  by  the 
grid.  This  average  strain  was  assumed  to  be  the  strain  at 
the  midpoint  of  the  grid.  The  distance  of  the  midpoint  from 
the  crack  tip  varied  from  0.041  to  0.294  inches  at  half  crack 
lengths  ranging  from  0.151  to  0.795  inches  for  the  aluminum 
and  from  0.195  and  1.00  inches  for  the  plexiglas. 
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EXPERIMENTAL  RESULTS 


Physical  observations  of  the  crack.  -  There  was  an  ob¬ 
jectionable  amount  of  surface  roughening  ahead  of  the  crack 
tip  for  the  longer  crack  lengths  in  the  6063-T6  aluminum 
shells.  This  roughening  made  it  difficult  to  find  the  exact 
location  of  the  crack  tip.  The  surface  roughening  was  worse 
in  the  6063-T6  aluminum  than  in  the  2024-T3  aluminum.  There 
was  a  negligible  amount  of  branching  of  the  crack  that  took 
place  during  its  growth.  What  branching  did  occur  quickly 
disappeared. 

Bulging.  -  Bulging  here  is  defined  as  large  outward  normal 
displacements  of  the  shell  in  the  crack  region.  Since  the 
elastic  stress  analysis  used  to  correlate  the  crack  propagation 
results  is  based  on  the  small  deflection  theory,  excessive 
bulging  may  reduce  the  reliability  of  the  final  results.  Hence, 
measurements  of  the  shell  surface  ahead  of  the  crack  in  shells 
#9,  #10  and  #11  were  taken  with  a  dial  indicator  to  see  how 
much,  if  any,  bulging  was  present.  Figures  9  and  10  are  plots 
of  the  radial  displacement  of  the  surface  of  the  shell  due  to 
the  presence  of  the  crack  only.  The  displacements  ahead  of  the 
crack,  caused  by  an  increase  in  applied  stress  of  6,900  psi, 
are  given  in  Figure  9  for  four  different  crack  lengths.  The 
curves  of  the  points  can  be  thought  of  as  profiles  of  the  shell 
surface  ahead  of  the  cracks.  Figure  10  compares  the  radial  dis¬ 
placement  at  the  crack  tip  of  shells  #9  and  #10  with  that  of 
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shell  #11  due  to  pressure  changes  of  13,500  psi  and  6,900  psi, 
respectively.  Since  the  maximum  operating  stresses  involved 
in  the  fatigue  tests  were  in  the  range  of  7,000  to  14,000  psi, 
these  curves  were  considered  as  the  extremes  of  the  radial 
di spl a  cements  for  all  specimens.  For  x<2.2,  the  radial  dis¬ 
placements  of  the  shell  surface  can  be  seen  to  be  less  than 
0.005  inches.  For  an  8  inch  diameter  shell,  this  may  be  con¬ 
sidered  as  "small  deflection"  and  hence  for  the  stress  range 
under  consideration,  bulging  does  not  appear  to  be  a  problem. 

Analysis  of  crack  growth  data.  -  The  proposed  crack  growth 
model  for  combined  loading  is 

2a-,  2(a-,+a9) 

da/dn  =  B(l+n)  (35) 

Substituting  from  (24)  and  (25)  into  (35),  we  obtain 
2a-,  2(a-,+a9) 

da/dn  =  B(l+fi)  [(Pg  +  yP  )K^g]  '  ^  (36) 

In  the  experimental  data  given  in  Appendix  II,  a°°,  2a  and 
An  are  known,  da/dn  is  obtained  by  differentiating  the  data 
by  means  of  a  three  point  central  difference  technique  (except 
at  the  ends,  where  a  forward  or  backward  difference  scheme  is 
used).  Using  (24),  (25),  (30)  and  (31),  the  stress  intensity 
factors  may  be  evaluated.  In  the  computer  program,  least 
squares  fit  of  log  da/dn  vs.  log  was  used  to  establish  the 
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values  of  B,  2ai  and  2(0-1+02)  for  y  “  1  sod  y  -  0.5.  The  re¬ 
sults  are  tabulated  in  Table  V  for  +  0.5  and 
for  the  range  of  convergence  of  the  theoretical  solution 
(X£2.2).  The  range  of  X  for  the  experimental  work  exceeded 
the  range  for  the  theoretical  solution.  However,  this  was  no 
great  handicap  since  the  growth  rate  for  x>2.2  was  very  rapid 
and  only  a  small  percentage  of  the  total  life  of  the  shell  was 
left.  The  term  r-j  in  Table  V  is  the  coefficient  of  correlation 
and  is  essentially  equal  for  +  0.5  K|^  and  +  K|^  for  each 
shell.  When  r-j  is  close  to  1.0,  as  it  is  for  the  shells  tested, 
the  assumption  of  a  linear  relationship  for  the  data  is  valid. 

The  term  S  is  the  standard  estimate  of  error  and  establishes 
y  X 

the  scatter  band  that  cover  95%  of  the  population. 

As  can  be  seen  in  the  table,  the  differences  in  A',  B,  Za-^ 
and  2(ai+a2),  for  twice  the  effect  of  bending,  are  too  small 
to  be  significant,  meaning  that  the  effect  of  bending  has  very 
little  influence  on  the  crack  growth  for  the  geometry  and  load¬ 
ing  conditions  considered.  For  the  remainder  of  the  analysis 
of  experimental  results,  y  =  0.5  is  used. 

jT  On  the  basis  of  combined  experimental  data,  which  is  given 
in  Figure  19,  the  crack  propagation  rate  for  6063-T6  aluminum 
may  be  expressed  as 


da/dn  =  1.068  x  10  (1+fi)  K 


(32) 


J  / 
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Figures  11  through  18  show  da/dn  vs,  =  K  +  0.5  K  . 

r  re  rb 

for  individual  shells  with  their  respective  scatter  bands. 

In  Figure  (19),  showing  the  combined  data,  the  three  basic 
levels  of  Q  are  designated  separately.  For  increasing  q,  the 
shift  of  the  data  to  the  left  is  characteristic  of  flat  plate 
fatigue  experiments. 

Figure  20  shows  the  comparison  of  the  crack  growth  rates 
for  Q  =  1.07  in  6063-T6,  2024-T3  and  7075-T6  aluminums.  The 
crack  growth  rates  for  2024-T3  and  7075-T6  aluminums  used  for 
the  figure  were  taken  from  [16]: 

da/dn  =  2.68  x  1  O' ^  ^  ( 1 +f2)  ^  K^)  ^  for  2024-T3 

da/dn  =  6.22  x  1 0'^°  ( 1 ^  ( K^)  ^  for  7075-T6 

The  line  for  6063-T6  aluminum  falls  between  those  of  2024-T3 
and  7075-T6  aluminums  for  growth  rates  greater  than  10'^ 
inches/cycle.  The  slope  of  the  line  is  almost  the  same  for 
6063-T6  and  7075-T6  aluminums. 

Transition  zone.  -  The  transition  zone  of  flat  to  shear 
fracture  gives  another  means  of  comparing  the  various  types 
of  aluminums.  According  to  Wilhem  [24],  the  zone  of  transi¬ 
tion  from  flat  to  shear  fracture  is  dependent  on  the  stress 
intensity  factor.  It  is  stated  in  [24]  that  for  2024-T3  and 
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7075-T6  clad  aluminums,  the  transition  takes  place  within  the 
range  of  of  3540  ±  1  770  #/inch^'^^.  To  see  if  this  range 
also  applies  to  6063-T6  aluminum,  the  range  stress  intensity 
factors  at  the  beginning  and  end  of  transition  of  the  fatigue 
cracks  of  five  specimens  were  determined  from  the  crack  lengths 
at  these  extremes.  For  the  purpose  of  calculating  the  range 
stress  intensity  factor,  it  was  assumed  that  the  shear  frac¬ 
ture  was  a  function  of  the  outer  surface  stress,  hence 

K  =  (Pu+P  )K  .  The  transition  initiation  and  termination 
r  '  b  e  re 

are  given  in  Table  IV.  The  results  for  the  6063-T6  aluminum 
fall  within  a  range  of  of  3440  ±  1220  #/inch^^^  which  is 
within  the  acceptable  range  of  2024-T3  and  7075-T6  aluminums. 

Strain  measurements  ahead  of  the  crack.  -  For  a  partial 
verification  of  theoretical  results,  static  strain  measure¬ 
ments  were  made  ahead  of  the  fatigue  crack  for  various  crack 
lengths  and  load  levels.  The  measured  strains  are  the  cir¬ 
cumferential  strains  in  an  internally  pressurized  cylinder 
with  fixed  ends,  which  contains  a  longitudinal  crack  free  of 
edge  loadings.  Thus  away  from  the  crack,  the  principal  strains 
are 


l-v2 

E 


00 

O 


(38) 


In  the  shell  loaded  at  the  crack  surface  only,  at  y  =  0  and 
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and  for  small  values  of  r  =  x 


•k 

a,  we  have 


0 


X 


(39) 


where  P  -  +  Pg*  Adding  circumferential  strains  e  given 

by  (38)  and  (39),  we  obtain  an  approximate  expression  for 
the  strain  ahead  of  the  crack  in  terms  of  stress  intensity 
ratio  P  as  follows: 


e 


ys 


c» 


(1  + 


p 

1+v 


(40  ) 


The  measured  strains  for  three  cylinders  (#12,  #14,  and 
#18)  as  well  as  that  obtained  from  equation  (40)  for  v  =  1/3 
are  shown  in  Figures  21  ,  22,  and  23.  Even  though  the  approxi¬ 
mate  strains  given  by  (40)  have  the  same  trend  as  the  measured 
strains,  the  agreement  is  obviously  unsatisfactory.  The  dif- 

k  k 

ference  is  due  to  finite  gage  length  ,  partial  plastic  defor¬ 
mations,  and,  perhaps  most  importantly,  the  fact  that  the  stress 
intensity  factor  (i.e.,  the  singular  term  alone)  does  not  cor¬ 
rectly  predict  the  perturbed  strains  caused  by  the  crack  in  a 


* 

Here,  for  the  results  of  a  higher  order  bending  theory 

(e.g.,  [22])  rather  than  equation  (22)  are  used.  In  the 
higher  order  theory,  which  satisfies  all  the  physical  boundary 
conditions,  ahead  of  the  crack  and  the  angular  dis¬ 

tribution  of  stresses  is  the  same  as  that  in  plane  extension. 

In  these  measurements,  1/16  inch  square  foil  gages  were  used. 
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homogeneous  strain  field.  In  the  shell,  the  details  of  the 
evaluation  of  stresses  away  from  the  crack  are  very  cumbersome 
and  difficult.  However,  a  reasonable  estimate  of  this  far 
field  effect  may  be  given  by  using  an  analogy  with  the  flat 
plate  solution. 

In  the  case  of  flat  plate  containing  a  crack  of  length 
2a  and  subjected  to  uniform  compressive  tractions  a“  on  the 
crack  surfaces,  the  strain  at  y  =  0  and  a  distance  r  from  the 
crack  tip  may  be  expressed  as. 


e 


yp 


''  (^==  - ') 
/r ( 2a+r ) 


(41  ) 


For  small  values  of  r,  T  1  and  (41)  gives  the  plate  equiva¬ 
lent  of  (39).  Thus,  if  we  assume  that  the  perturbed  stress 
fields  in  plate  and  shell  are  similar  up  to  distances  r  under 
consideration,  an  approximate  value  of  for  the  shell  may 
be  obtained  by  multiplying  (39)  by  factor  T  which  is  given  in 
(41).  With  this,  (40)  may  be  modified  as 


E 


ys 


00 


(1  + 


P*T 

1+v 


(42) 
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Again,  for  v  =  1/3,  the  comparison  of  measured  strains  with 
those  obtained  from  (42)  is  seen  in  Figures  24,  25  and  26. 

Since  T  is  always  less  than  unity,  (40)  gives  a  higher  (theo¬ 
retical)  estimate  for  the  strains  -  partly  explaining  the  dis¬ 
crepancy  in  Figures  21  to  23.  However,  differences  also  exist 
in  Figures  24  to  26.  This,  of  course,  is  due  partly  to  the 
assumption  concerning  the  similarity  of  stress  fields  in  flat 
plate  and  shell  which  leads  to  (42),  and  partly  to  inelastic 
deformations  around  the  crack  tip.  The  fact  that  the  differ¬ 
ences  are  more  noticeable  at  higher  load  levels  (i.e., 
and  larger  a/r  ratios  substantiates  the  validity  of  the  latter 
content! on . 

The  effects  of  a  finite  gage  width  can  be  estimated  by  ex¬ 
amining  the  error  involved  in  assuming  the  mean  strain  to  be 
at  the  center  of  the  gage.  As  seen  from  Figure  27,  the  actual 
mean  strain  is  closer  to  the  crack  tip  than  the  center  of  the 
gage.  The  mean  strain  point  is  established  when  the  shaded 
areas  are  egual.  This  means  that  the  data  will  be  plotted 
closer  to  the  ordinate  than  it  should  be.  This  puts  the  data 
above  the  prediction  curve  for  the  higher  applied  strain  or 
smaller  distances  r.  For  small  applied  strains  or  large  dis¬ 
tances  r  where  the  elastic  strain  curve  is  relatively  flat, 
the  measured  strain  and  the  mean  strain  are  essentially  equal, 
and  the  effect  of  the  finite  gage  length  is  negligible  (see 
Figure  24). 
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To  obtain  further  knowledge  of  the  stresses  in  the  crack 
tip  area,  static  tests  were  conducted  on  Plexiglas  G  of  a  plate 
under  tension  and  a  cylinder  under  pressure.  Strains  were 
measured  ahead  of  the  crack  at  room  temperature.  Plexiglas  was 
chosen  so  that  the  distance  of  a  strain  gage  ahead  of  the  crack 
tip  could  be  measured  when  the  gage  was  mounted  on  the  inner 
surface  of  the  cylinder.  The  nominal  dimensions  were  12"  x 
24"  X  1/8"  for  the  plate  and  9  1/2"  diameter,  20"  length  and 
1/8"  thickness  for  the  cylinder. 

The  pressure  in  the  cylinder  was  measured  with  a  dead 
weight  tester  since  the  test  pressures  were  too  low  to  be  ac¬ 
curately  controlled  with  the  test  set  up  described  in  this  re¬ 
port.  Higher  test  pressures  were  not  possible  without  risking 
catastrophic  failure. 

The  stress  intensities  on  the  inner  and  outer  surfaces  of 
the  shell  were  computed  from  the  strain  readings  ahead  of  the 
crack.  They  were  then  compared  with  the  stress  intensity  com¬ 
puted  from  the  strain  reading  ahead  of  an  equal  length  crack 
and  applied  stress  in  a  plate.  The  comparison  was  made  for 
three  different  crack  lengths.  The  ratio  of  the  stress  inten¬ 
sities  in  the  cylinder  to  that  of  the  plate  should  equal  P^  + 

Pb  Pg  "  outer  and  inner  surface  of  the  cylinder, 

respecti vely . 
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Figure  28  is  a  plot  of  the  ratio  of  the  stress  intensity 
factors  vs.  crack  length  for  the  test  on  the  plexiglas  and  for 
the  eight  terms  and  Folias'  [2]  solutions.  These  experimental 
results  are  largely  qualitative.  They  were  not  pursued  further 
because  of  the  steady  drift  in  the  gage  readings  caused  by 
creep  and  the  fact  that  the  creep  compliance  for  the  two  mate¬ 
rials  were  not  quite  the  same. 
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CONCLUSIONS 


I  The  results  obtained  in  this  study  indicate  that  within 
the  practical  range  of  fatigue  crack  propagation  rates,  i.e., 

10  to  10  in/in/cycle,  the  stress  intensity  factor  is  an 
effective  correlation  parameter  in  analyzing  the  fatigue  data, 
and  provides  a  powerful  means  of  predicting  fatigue  crack 
growth  characteristics  of  shells  from  that  of  flat  plates. 

With  the  exception  of  problems  involving  "low  cycle  fatigue", 
from  the  design  engineer's  view  point,  essentially  the  problem 
may  be  reduced  to  one  of  evaluating  the  stress  intensity  fac¬ 
tor.  After  this,  one  may  use  the  flat  plate  results  correlated 
against  the  stress  intensity  factor  to  predict  the  proper 
growth  rates. 


For  the  "low  cycle  fatigue",  it  is  suggested  that  a  crack 
propagation  model  based  on  the  plastic  deformations  around  the 
crack  tip  mentioned  in  this  report  may  prove  to  be  more  effec- 
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APPENDIX  I 


Details  Of  Elastic  Solution 
The  coupled  differential  equations  to  be  solved  are 


Eha^  32W 
R 


9X 


+  V^F  =  0 


v‘*W 


a2  3  2f 

RD 


0 


(I.l) 


(1.2) 


subject  to  the  following  boundary  conditions  on  the  crack  sur¬ 
face 


My(x,0)  = 


D  r92W  .  a^W-i  _  '  0 


Dm 


\i  f  n\  -  r  9  j.  f  o  'i  9  1  _  0 

Vy(x,0)  =  -  (2  -  v)  gx^ 9y ^  ”  a ^ 


N,y(x,0)  = 


1  a^F  ^  \ 
"a^  axay  a 


(1.3) 

(1.4) 

(1.5) 

(1.6) 


and  the  continuity  requirements 


lim  (-^  W'^  -  W)  =  0,  n  =  0,1, 2, 3 
|y|->o  9y  9y 


+  a" 

lim  (^-  F^  -  ^  F  )  =  0,  n  =  0,1, 2, 3 
|yho  9y"  ay 


(1.7) 


(1.8) 
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along  the  prolongation  of  the  crack.  Also  W  and  F  and  their 
first  derivatives  must  be  finite  far  away  from  the  crack. 

For  the  symmetric  case,  i.e.,  v  =0  and  t  =  0,  W  and  F 

0  0 

can  be  represented  as 


00 


W(x,y) 

=  /  w(s,y)  cos  xs  ds 

0 

(1.9) 

F(x,y) 

00 

=  /  f(s,y)  cos  xs  ds 

(1. 10) 

0 


If  uniform  convergence  is  assumed,  then  the  order  of  integration 
and  differentiation  can  be  interchanged.  Therefore 


v'^F  =  /  [s'*f  -  2s2  0  +  0]  cos  xs  ds 

v'^W  =  /  [s‘»w  -  2s2  cos  xs  ds 

=  j  ~  s^f  COS  xs  ds 
0 

7  2  ^ 

— =  /  -  s^w  cos  xs  ds 
0 

Equations  (I.l)  and  (1.2)  can  now  be  written  as 


/  [-  s2w  +  s^^f  -  2s2  ^  cos  xs  ds  =  0  (1. 11) 

/  [s'»w  -  2s2  0  +  0  +  f]  cos  xs  ds  =  0  (1.12) 
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For  the  integrals  to  be  zero,  the  integrands  must  vanish.  There¬ 
fore,  we  have 


R  i 


d'+w  0-9  d^w  .  -u^  ,  f  _  o 


the  solution  of  which  may  be  written  as 


8  r  .y 
w  =  ^  Q.e  J 

1  J 


8  r  .y 
f  =  I  R^e  J 
1  J 


where  Q.  and  R*.  are  arbitrary  constants  and  r.  are 
J  J 


r-j  =  +  /s(s  -  aX) 
T2  -  -  >^s(s  -  aXj 
Fj  =  +  /S.(S-  +  aX) 

=  -  /s(s  +  aX) 


(1.13) 

(1.14) 

(1.15) 
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Fg  =  +  /s(s  -  ex) 


Fg  =  -  /s(s  -  3  A) 

(1.16) 

Fy  =  +  /s(s  +  ex) 

Fg  =  -  /s(s  +  ex) 


„  =  /T.  B  =  /IT,  x** 


12(l-v2)a‘' 


For  |y|-^“,  w  and  f  must  vanish;  therefore,  Q.  =  R'.  =  0  for  odd 

J  J 

j.  The  remaining  constants  are  not  linearly  independent.  Thus, 
substituting  equations  (1.15)  into,  say  (1.13),  we  obtain 


r;.  =  [ 


- ! - Z-]  Q,- 

(s2  -  r?)  J 


(1.17) 


Now  substituting  the  values  for  f.  from  equations  (1.16)  into 

J 

(1.17),  we  find 


=  -  i  /EhD  Q2^4 
R^,8  =  i  /EhD  Qg^g 


(1.18) 


Letting  Q2,4,6,8  ^1,2, 3,4 


from  equations  (1.9), 


(I.IO), 
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(1.15)  and  (1.18),  we  obtain 


W(x,y-)  =  /  (P^e"*^^^^  ■  1^1 

+  p  g-/s(s  -  BA)  |y| 

3 

+  p^e-'^s(s  +  BA)  |y|j  (1.19) 


and 


F(x,y")  =  -  1  /EhD  /  (p  lyl 

0 

^  P  g-/s(s  +  aX)  |y|  _  p  g-/s(s  -  BX)  |y| 

2  3 

_  p^g-/iTs“TTrr  |y|)  ggg  (1.20) 

where  ±  sign  refers  to  y>0  and  y<0,  respectively.  Equations 
(1.19)  and  (1.20)  are  the  same  as  those  obtained  by  Folias  [2]. 
The  procedure  to  be  followed  to  arrive  at  the  corresponding 
singular  integral  equations  is  then  given  in  [2].  Thus  defining 


00 

Ui(x)  =  /  [  /s  ( S  -  aX )  P"!  +  /S  ( S  +  aX )  ?2 
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+  (/s(s  -  aX)  P-J  -  /s(s  +  aX}  P2)]  COS  xs  ds 


(1.21) 


U2(x)  =  /  [/s(s  -  aX)  P-j  +  /s  (S  +  aA)  P2 


^  (/s(s  -  Oix)  P.| 


-  /s(s  +  aX)  P2)]  COS  XS  ds 


X  <1  ,  y  =  0  ) 


(1.22) 


we  obtain 


2iTn 


/  [u,(?)L,  +  =  - =  x;  |x|<l,  |5|<1  (1.23) 

1  '  '  ^  ^  i  /EhD 


/  [u^(?)L3  +  U2(5)L4]  d§  =  -  |x|<1,  U|<1  (1.24) 


where  the  kernels  L-j  ,  L2»  L^,  and  are  (c  =  x  -  c) 


o(  ^  C  ^ !/  (  3  ^  I  C 


L,  =  aX  s1n  ( 


)  +  cos 


^6X  s1n  (^)K;,  (2i^)  t  cos  (^)K,  (^4^) 


+  V 


2a  C 


sin  (2^)K,  (^ii) 


TsTn  (  2  ^“^l  ^  2  ^ 


(1.25) 


45 


L2  =  (v^  -  1)  BX  sin  (^)K^  (^4^) 

+  (vq  -  1)  4f  ''OS  (^)Kl  (^4^^ 

"  2^  '*"  <^>^1 

"2bU|  (4^)'<T  (^4^)  (^*26) 

L3  =  (v^  -  |)  ax  [sin  (4^)K^  (^^)] 

+  (vq  -  |)  4f  tcos  (^)K^  (^4^)^ 

+  6X  (-  v2  +  Vq  -  j)  [sin  (^)Kq  (—4  ^  )] 

ftf  ^0  ^  (^4^)] 

-  ^0  ^  tsin  (^)K^  (iA^)] 

^  =  fftr  (4^)'<i  (^4^)] 

+  fA  (4^^  _  3  .  2v2)  [sin  (^)K^  (^^)] 
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+  fjfj  ('I'-o  -  H  -  I’ 

*  2?W  -  "o)  [  =  '"  <¥>^1 

^  2^  [  =  '^>^1 

-  [s1n  (2ii)K^  (^4^)]  (1-28) 

In  the  foregoing  equations,  ^  ~  ^  *^n  ^  Bessel 

function  of  the  third  kind  of  order  n. 

The  solutions  u^(x)  and  U2(x)  must  be  Holder  continuous 
for  some  positive  Holder  indices  y-i  ,  and  y2  ^  ^ 

closed  interval  (-1,1),  and  Ui(x)  and  U2(x)  are  bounded  near 
the  ends  of  the  crack.  These  details  are  discussed  in  [1]. 

The  kernels  can  now  be  written  in  series  form.  Folias 
eliminates  all  terms  of  order  and  greater  in  his  series  ex¬ 
pansions  which  puts  a  severe  limitation  on  the  useful  range  of 
X.  In  this  report,  polynomial  approximations  of  the  Bessel 
functions  are  used  which  contain  terms  up  to  X^*^,  and  are  good 

for  x<2.  The  sine  and  cosine  series  are  expanded  up  to  terms 

-  8  “9 

of  x^^  with  corresponding  errors  of  2.5  x  10  and  3.2  x  10 
respectively,  for  x<2.  This  is  in  the  range  of  error  for  the 
Bessel  functions.  The  Bessel  functions  are  regular  throughout 
the  (x-5)  plane,  cut  along  the  negative  real  axis.  Multiplying 
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all  series  together  that  are  now  present  In  the  kernel  and 
collecting  terms,  the  integral  equations  can  now  be  written  in 
the  form 


i2TrnQX  1  (1  +  Vq)u^(5)  +  (vq-1)u2(c) 

(EhD)‘/a= 

1 

+  X  /  {[(1+v  )Ui(c) 

-1  °  ' 

+  I  H  x'">*3,,.5)4b*3 

+  i[u,(5)  +  UpCc)]  [(1-vJ  [p 

*  ^  u  Q  in 

+  I  log 

*  '’o  [f 

+  I  log  ^1^]]}  ds  (1.29) 
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and 
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(1.30) 


.  I  log  de 

0  D 

where  the  constants  D^,  etc.,  are  listed  in  Appendix  III. 
These  constants  are  independent  of  the  shell  parameters  and 
materials  constants. 

The  unique  solution  of  (1.29)  and  (1.30)  are  of  the  form 


°°  O 

u,(?)  =  /r^  I  A  1?I<1 


(1.31) 


p=0 


Uj,(?)  =  I  B  x^P(l-c2)P;  U|<1 

p=o  ^ 


(1.32) 


where  Ap  and  Bp  are,  generally,  complex  functions  of  A. 

Substitution  of  (1.31)  and  (1.32)  into  (1.29)  and  (1.30) 
leads  to  the  following  integrals  which  have  to  be  evaluated: 


I  (x)  = 
P 


C.P.V. 


I  X-E 


log  |x-c|  d? 

Pn  _  ] 


(1.33) 


where  n  is  an  odd  integer. 
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Through  contour  integration  and  using  a  table  of  integrals, 


we  obtain 


Ip(x) 


TT  Y  (  -  1  ) 

s  =  o  2!2^ 


n  (2p+3-2b)x^P^^"^^ 
b  =  o 


(1.34) 


(x) 


TT 


P 

n 

s  =  o 


(2p+l-2s) 

(2p+2-2s) 


r=  1 


n ! 

(n-2r) ! (2r)  ! 


n 

b  =  l 


(2b-l)  ..n 
(2p+2+2b)^ 


2r 


(1.35) 


where  q  =  (n-l)/2  for  odd  n  and  q  =  n/2  for  even  n.  '^pn(^) 
can  be  determined  by  first  considering,  for  n  =  0, 


'^00^^)  "  ^  (1-?2)P+''/2 

pu  _  1 

Differentiating  with  respect  to  x,  we  have 
=  C.P.V.  I  (1-§2)P+1/2 

where  Ip(x)  is  known.  Integrating  Ip(x)  in  x  gives  JpQ(x)  plus 
a  constant  of  integration,  or 
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Jpo'x'  = 


TT 


2pt3 


s  =  o 


(-1)^ 

s!2^ 


s 

n 

b  =  o 


(2p+3-2b)  2p+2-2s 
{2p+2-2s)^ 


(1.36) 


Jnn(O)  =  /  log  ui  dc 

P*^  -  1 


p!(-1)'^(2r+n-1)! 


b  P!  V"  I  ;  \ci-ru-  I  }  : _ r  ■ 

'  Jo  r!(p-r)!2''*"/2*l(,,„/2*i)! 


2r+n  /  ,\k 

+  log  2  +  I  — — ] 

k=l 


for  n  even 


(1.37) 


Jpn(O)  =  0,  for  n  odd. 

Letting  now  n  =  1  in  Jpp(x)  aod  differentiating  with  re¬ 
spect  to  X,  we  find 


J'(x)  =  /  +  j  log  |x-{|  dE 

pi  -1  -1 


But  from  (1.35) 


P  (2p+1-2s) 
^  (2p+2-2s) 


Hence , 

JMx)  .  Ip„(x)  +  Jp„(x) 
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Or,  in  general,  we  obtain 


'^pn<’‘)  =  !  ‘p(n-l)<’’)‘'>‘  Jp(n-1)<'‘>‘'’' 


+ 


Jpn<»> 


(1.38) 


Therefore ,  J 
Ip„(x)  and  J 


pn 

pn 


(x)  can  be  reduced  to  a  combination  of  I  (x), 

p 

(0). 


With  the  i ntegrals  gi ven  by  Q.33)  known,  (1.29)  and  (1,30) 
can  be  written  in  series  form  in  odd  powers  of  x.  The  contri¬ 
butions  of  I  (x),  I  (x)  and  Jp,_(x)  to  each  coefficient  of 
2  1 

x^  ,  where  m  =  0,1,...,  13,  are  given  below. 


For  I  (x) 
P  ' 


X 


TT 


(-l)P 

p!2P 


P 

n 

s  =  o 


(2p+3-2s) 

2p+3 


X3  IT 


(-l)P'^  P;^  (2p^^3-2s) 

(p-l)!2P“^  s  =  o 


2m+l  .  (-l)P'"’  P'"’  (2p+3-2s) 

A  TT  - - — —  n  k To - 

(p-m)!2P'"'  s  =  o  ^P"^^ 
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X  TT 


P 

n 

s  =  o 


(2p+l-2s) 

(2p+2-2s) 


+ 


C^(2n+1) 


n 

n 


b  =  l 


{2n+l-2b)i 

(2p+2+2b)-' 


^  IT 


P 

n 

s  =  o 


(2p+l-2s) 

(2p+2-2s) 


13 


(2n+l)(2n)(2n-l)  "  (2n+l-2b) 

3!  (2p+2b) 


X®  tt 


P 

n 

s  =  o 


(2P+1-2S) 

(2p+2-2s) 


13 

[Co  +  I  C 
^  n=3 


(2n+l)'’'(2n-3) 

5! 


n 

n 

b  =  3 


(2n+l-2b)n 

(2p-2+2b)-' 


x2«tl 


H 


P 

n 

s  =  o 


(2p+l-2s) 

(2p+2-2s) 


[C 


m 


13 

I , 

n=ni=l 


2ni+l 

n  (2n+r-2m) 


(2m+l)!  b=m+l 


(2n+l-2b)  n 
(2p+2+2b-2m)-' 


where  are  appropriate  constants  E^,  etc. 


For  J  (x) 
pn 


X  - 


y  C  P'  (2s  +  2n-1)!  !(2n+n 

L  n  cfT*rrZ”c~TT  r*  ir»  _L  1 


s  =  o  n  =  o  "  S!(p-s;!  2^''"^^(s+n+l)!  2s  +  2n+2 


2s+2n  /  nxk 

log  2  +  I  }  + 


P 

n  n 


2p+2-2s  ‘-^0 


1 3  n 


13 

letting  \ 
n  =  1 


^n^(n-l)p  =  hp 


x3  -v  B,  +  TT  ^  n  ^  )  _jn P  (2p+l-2s)  [-op 

p!2P  s  =  o  °  (2p+3)2.3  ^  2.3  {2p+2-2s)  l-^^l 


I  c  (20+1)  ; 

n=2  ^  b=2  2p+2b 


+  TL  P  <:2p+1-2s)  rr  +  r  (2n)(2n-l)  JJ  (2n+l-2b)-, 
3  l2p+2-2s)  Lh  +  ^ --(2p+2b)  ^ 


y  C„D  +  C  E 
n=i  1  np  0  np 
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ii 


^n'^(n-l  )p  "  ^2p 


13 


2 


2p 


C  ^(-1)P~1 . 

°  (p-l)!2'^“^  s  =  o 


(2p+3-2s)  , 
(2p+3)4.5 


(-I)P  P  (2p+3-2s)2.3 
p,2P  s=o  ^2p.3)5! 


+ 


IT 

4.5 


P 

n 

s  =  o 


(2p+l-2s) 

(2p+2-2s) 


[5C 


2 


.  'I  C„(2n*l) 
n=  3 


n 


n 

b=3 


(2n+l-2b)n  ^  ^ 
(2p+2b-2)-'  5 


P 

n 

s  =  o 


(2p+l-2s) 

(2p+2-2s) 


[C 


2 


.  „  (2n)(2n-1)(2n-2)(2n-3)  "  ( 2n-»-1 -2b ) -, 

n=3  "  b  =  3  (2p+2b-2)-' 


13  1 

y  c  F  +  y  c  G 

n^2  "  "P  n=o  "  "P 
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etc.,  where  are  appropriate  constants  H|^,  L|^,  etc. 

By  comparing  coefficients  of  x  in  the  resultant  series  form 
of  (1.29)  and  (1.30),  a  system  of  2p  +  2  simultaneous  equations 
can  be  formed  to  solve  for  Ap  and  Bp.  As  the  number  of  simul¬ 
taneous  equations  increases,  the  value  for  a  particular  Ap  or 
Bp  converges  to  its  exact  value.  Since  the  stresses  at  the  crack 
tip  are  desired,  it  is  necessary  only  to  compute  the  coefficients 

A„  and  B^. 

0  0 

The  computer  programs  for  comparing  coefficients  are  given 
in  Appendix  IV.  For  the  first  program  (coefficients  of  x”)  the 
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range  and  increment  of  A  must  be  established  along  with  the 
maximum  number  of  simultaneous  equations  desired.  The  card  out¬ 
put  of  the  first  program  becomes  the  input  of  the  second  program. 
For  the  second  program  (solution  for  and  P||^)  the  values  of 
V  and  m^  must  first  be  established.  This  way,  the  output  of  the 
first  program  can  be  used  for  various  values  of  v.  The  output 
of  the  second  program  gives  the  values  of  A^,  Pjl^  and  P^. 

P^  and  Pg  are  designated  as  PlOA  and  RP20A,  respectively,  in 
the  program. 

From  (21-24)  and 

oo 

/  s"^J  (gs)cos  xs  ds  = 

0 

/^(2g)-^^[r(w+l/2)]-‘'(g2 

{ 

0 


-x2)^-''/2.  o<x<g 

;  Rep>-  1/2} 

;  g<x<«> 
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/s (s-xe)  P3(s) 


V  S  1 

/  0  I  \  0  0  I 

^3X  ■  2^  2  s 


A„-B^  J,(s) 

g  /  0  0  \  1  . 

23  ^  2  ^  gX 


(1.41) 


/s(s+X3)  P4(s) 


V^s  ,  A„+B„  Ji(s) 

■  (if-  *  1)  (-V^>  -V- 


A„-B^  Ji(s) 

a _  /  0  0  \  1 

23  '  2  '  gX 


(1.42) 


The  stresses  may  then  be  obtained  by  a)  substituting  equations 
8:39)  through  (142)  into  equations  (19)  and  (20);  b)  substituting 
the  result  into  the  boundary  condi ti ons, equati ons  (3)  through 
(6);  and  c)  changing  the  order  of  integration  and  differentia- 
ti  on . 
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APPENDIX  II 


crack  Growth  Data 

Shell  #9  -  6063-T6  Shell  #10  -  6063-T6 


Mean  stress  =  6850  psi  Mean  stress  =•  7220  psi 

Stress  range  -  6300  psi  Stress  range  ■=•  6770  psi 


2a 

An 

2a 

An 

0.28620 

0,000 

0.31280 

0,000 

0.30610 

3,000 

0.34320 

3,000 

0.32310 

2,000 

0.36270 

2,000 

0.33470 

2,000 

0.38400 

2,000 

0.35820 

2,000 

0.41310 

2,000 

0.38300 

2,000 

0.44970 

2,000 

0.40640 

1,500 

0.49010 

1,000 

0.42420 

1,000 

0.52250 

1,000 

0.43060 

1,000 

0.55050 

1,000 

0.44905 

1,000 

0.60010 

1,000 

0.46365 

1,000 

0.63720 

600 

0.47835 

1,000 

0.69000 

600 

0.49540 

1,000 

0.75670 

600 

0.51085 

1,000 

0.82720 

600 

0.53320 

1,000 

0.89240 

300 

0.55655 

1,000 

0.94890 

200 

0.58300 

1,000 

1.03220 

200 

0.62025 

1,000 

1.11260 

100 

0.64550 

600 

1.17960 

100 

0.67170 

600 

1.21130 

100 

0.69660 

600 

0.74250 

600 

0.76110 

300 

0.79300 

300 

0.83430 

300 

0.87125 

300 

0.91830 

300 

0.96850 

300 

1.00325 

200 

1.05730 

200 

1.08470 

100 

1.12240 

100 

1. 15820 

100 

1.20280 

100 

1.26430 

100 

1.32980 

75 

1.38800 

45 

1.42360 

30 

1.46950 

30 

1.49250 

30 

1.56720 

30 

1.64260 

30 
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Shell  #11  - 

6063-T6 

Shell  #11 

cent. 

Mean  stress 

s  3590  psi 

2a 

An 

Stress  range  =  3370  psi 

0.51730 

6,000 

0.52220 

6,000 

2a 

A  n 

0.54410 

6,000 

0.28900 

0,000 

0.55770 

6,000 

0.29320 

6,000 

0.56610 

6 , 000 

0.29850 

6,000 

0.57430 

6,000 

0.29960 

6,000 

0.58610 

6,000 

0.30280 

6,000 

0.59560 

6,000 

0.30560 

6,000 

0.60380 

6,000 

0.30800 

6,000 

0.62840 

6,000 

0.31280 

6,000 

0.64720 

6,000 

0.31620 

6,000 

0.66160 

6,000 

0.32330 

6,000 

0.68450 

6,000 

0.32530 

6,000 

0.70430 

6,000 

0.32610 

6,000 

0.72880 

6,000 

0.33040 

6,000 

0.75840 

6,000 

0.33700 

6,000 

0.83490 

6,000 

0.34010 

6,000 

0.87040 

6,000 

0.34440 

6,000 

0.92110 

6,000 

0.34630 

12,000 

0.94750 

3,000 

0.34910 

6,000 

0.97910 

3,000 

0.35320 

6,000 

1.02800 

3,000 

0.35370 

12,000 

1.05350 

3,000 

0.36630 

6,000 

1.10690 

3,000 

0.36980 

6,000 

1.15780 

3,000 

0.37010 

6,000 

1.22710 

3,000 

0.37450 

6,000 

1.25570 

1,000 

0.37600 

6,000 

1.28530 

1,000 

0.38270 

6,000 

1.32240 

1,000 

0.38400 

6,000 

1.35580 

1,000 

0.38810 

6,000 

1.39430 

1,000 

0.39470 

6,000 

1.44950 

1,000 

0.40090 

6,000 

1.53580 

1,000 

0.40440 

6,000 

1.58210 

600 

0.40860 

6,000 

1.65430 

600 

0.41380 

6,000 

1.75300 

600 

0.42080 

6,000 

1.88270 

600 

0.42720 

6,000 

2.00240 

312 

0.43770 

6,000 

0.45300 

12,000 

0.46310 

6,000 

0.47300 

6,000 

0.47910 

6,000 

0.48430 

6,000 

0.49510 

6,000 

0.50460 

6,000 
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Shell  #12  - 

6063—T6 

Shell  #12 

cent. 

Mean  stress 

=  5370  psi 

2a 

/in 

Stress  range  =  5040  psi 

0.56040 

1,500 

0.56770 

1,500 

2a 

0.57930 

1,500 

0.33450 

0 , 000 

0.59360 

1,500 

0.34220 

1,500 

0.60710 

1,500 

0.34420 

1,500 

0.62290 

1,500 

0.35240 

1,500 

0.65300 

1,500 

0.35480 

1,500 

0.65980 

1,000 

0.36340 

1,500 

0.66890 

1,000 

0.37090 

1,500 

0.68900 

1,000 

0.37800 

1,500 

0.69620 

1,000 

0.33740 

1,500 

0.70560 

1,000 

0.39510 

1,500 

0.72300 

1,000 

0.40070 

1,500 

0.78370 

1,000 

0.41290 

3,000 

0.79750 

1,000 

0.41620 

1,500 

0.81930 

1,000 

0.41750 

1,500 

0.85630 

1,000 

0.42060 

1,500 

0.88090 

1,000 

0.42300 

1,500 

0.91360 

1,000 

0.42520 

1,500 

0.94950 

1,000 

0.42560 

1,500 

1.00120 

1,000 

0.42850 

1,500 

1.09680 

1,000 

0.42930 

1,500 

1.18510 

600 

0.42980 

1,500 

1.21030 

150 

0.43320 

3,000 

1.24630 

150 

0.43690 

1,500 

1.30530 

150 

0.44070 

1,500 

1.38350 

150 

0.44270 

1,500 

1.44970 

100 

0.44640 

1,500 

1.48700 

100 

0.44710 

1,500 

1.55740 

100 

0.44860 

1,500 

1.67910 

100 

0.45370 

1,500 

0.45730 

1,500 

0.45940 

1,500 

0.46560 

1,500 

0.46940 

1,500 

0.47580 

1,500 

0.47880 

1,500 

0.48920 

1,500 

0.49790 

1,500 

0.50750 

1,500 

0.51480 

1,500 

0.52130 

1,500 

0.53130 

1,500 

0.53970 

1,500 

0.55170 

1,500 
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Shell  #13  - 

6063-T6 

Shell  #14  - 

6063-T6 

Mean  stress 

=  11080  psi 

Mean  stress 

=  5460  psi 

Stress  range 

=  3400  psi 

Stress  range  •=  5070  psi 

2a 

^^n 

2a 

An 

0.51790 

0,000 

0.28800 

0,000 

0.52680 

2,000 

0.30380 

3,000 

0.53760 

2,000 

0.31530 

3,000 

0.54700 

2,000 

0.32990 

3,000 

0.55740 

2,000 

0.34250 

3,000 

0.58920 

2,000 

0.35640 

3,000 

0.59870 

2,000 

0.37490 

3,000 

0.60770 

2,000 

0.39010 

3,000 

0.61640 

2,000 

0.43630 

3,000 

0.62810 

2,000 

0.46130 

3,000 

0.64240 

2,000 

0.48830 

3,000 

0.65850 

2,000 

0.51920 

3,000 

0.67170 

2,000 

0.55080 

3,000 

0.68820 

2,000 

0.59870 

3,000 

0.70370 

2,000 

0.66260 

3,000 

0.72050 

2,000 

0.74190 

1,500 

0.73160 

2,000 

0.80430 

1,500 

0.74890 

2,000 

0.84530 

1,000 

0.76630 

2,000 

0.89250 

1,000 

0.78440 

2,000 

0.95310 

1,000 

0.81170 

2,000 

1.02590 

900 

0.83320 

2,000 

1.11490 

600 

0.85840 

2,000 

1.21180 

600 

0.88040 

1,500 

1.39390 

600 

0.90100 

1,500 

1.43850 

100 

0.92500 

1,500 

1.51110 

100 

0.95730 

1,500 

1.61260 

60 

0.98810 

1,500 

1.02400 

1,500 

1.05530 

1,000 

1.08860 

1,000 

Shell  #15  - 

6063-T6 

1.13280 

1,000 

1.17450 

1,000 

Mean  stress 

=8420  psi 

1.24830 

1,000 

Stress  range  =  5120  psi 

1.30170 

600 

1.40980 

300 

2a 

A  n 

1.44390 

300 

0.28130 

0,000 

1.48710 

300 

0.29500 

3,000 

1.53030 

300 

0.30880 

3,000 

1.56740 

200 

0.32070 

3,000 

1.60530 

200 

0.33390 

3,000 

1.62170 

100 

0.35100 

3,000 

1.64340 

100 

0.36650 

3 , 000 

1.67250 

100 

0.39070 

3,000 

1.69460 

100 

0.41720 

3,000 
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shell  #15  cont. 

Shell  #13 

cont. 

2a 

An 

2  a 

0.44460 

3,000 

0.56170 

3,000 

0.47680 

3,000 

0.58920 

3,000 

0.51440 

3,000 

0.62890 

3,000 

0.55910 

3,000 

0.67040 

3,000 

0.61560 

3,000 

0.72020 

3,000 

0.64430 

1,500 

0.75680 

1,500 

0.68910 

1,500 

0.80430 

1,500 

0.72560 

1,500 

0.85100 

1,500 

0.79160 

1,500 

0.93120 

1,500 

0.85290 

1,500 

1.02170 

1,000 

0.91370 

1,000 

1.17970 

1,000 

0.97620 

1,000 

1.37290 

600 

1.08060 

1,000 

1.59050 

150 

1.15850 

600 

1.24910 

600 

1.33430 

300 

1.42010 

200 

Shell  #19 

-  SAE  1018  Steel 

1.45870 

100 

1.48930 

60 

Mean  stress  =  7000  psi 

1.50910 

60 

Stress  range  •=  7000  psi 

2a 

An 

0.32575 

0,000 

Shell  #18  - 

6063-T6 

0.32635 

6,000 

0.32645 

6,000 

Mean  stress 

=  5350  psi 

0.32720 

6,000 

Stress  range  5070  psi 

0.32755 

12,000 

0.32800 

6,000 

2a 

0.32830 

6,000 

0.23820 

0,000 

0.33080 

6,000 

0.30200 

3,000 

0.33455 

12,000 

0.30970 

3,000 

0.33655 

6,000 

0.32320 

3,000 

0.33925 

4,000 

0.33930 

3,000 

0.35550 

3,000 

Mean  stress  =  7500  psi 

0.37220 

3,000 

Stress  ranqe=  7500  osi 

0.38340 

3,000 

0.39370 

3,000 

2a 

yy  n 

0.40900 

3,000 

6.33925 

00,000 

0.42080 

3,000 

0.34070 

12,000 

0.43850 

3,000 

0.34395 

12,000 

0.45520 

3,000 

0.34425 

12,000 

0.47400 

3,000 

0.34660 

12,000 

0.49440 

3,000 

0.34915 

6,000 

0.51450 

3,000 

0.35380 

6,000 

0.53510 

3,000 
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APPENDIX  III 


H  "  Vn4Wn-Yn 


-5.5791176x10”^ 

7.1585753x10“2 
■1. 5027655x10-3 
•6.1278490x10-3 
8.8787107x10-7 
6.3429898x10-5 
•6.0045804x10"^^ 

■1.6814260x10-^3 

1.1431171x10-^3 

1.5924912x10-^® 

■7.4005562x10-31 

■5.2590588x10-34 

1.2557491x10-36 

2.8131237x10-30 


-5.1318632x10-3 

7.1585759x10-3 

-1.4570399x10-4 

-3.4002697x10-6 

4.4362985x10-® 

2.4395960x10-10 

■2.1138323x10-13 

■4.9414796x10-13 

2.8896248x10-17 

3.6827901x10-30 

■1.0961288x10-32 

■9.0371123x10-26 

0.0 


-3.9269908x10-1 

1.1866800x10-3 

1.6106839x10-3 

-2.9019260x10-3 

-4.0143251x10-7 

-1.8654274x10-10 

1.9537601x10-11 

-1.5540472x10-13 

-3.0005865x10-16 

1.8927971x10-1® 

1.6359083x10-31 

-7.4051638x10-34 

-2.2267587x10-37 

0.0 


0 

1.9634954x10-1 

-9.8174770x10-3 

1 

-2.3651856x10-3 

-1.7050709x10-3 

2 

1.3422335x10-3 

2.1092260x10-4 

3 

-2.4743023x10-3 

-3.4679912x10-6 

4 

-3.6127755x10-7 

-3.4668622x10-® 

§ 

4.2211942x10-® 

3.6464548x10-10 

6 

1.8149490x10-11 

1.2562638x10-13 

7 

-1.4354576x10-13 

-9.0591396x10-13 

8 

-2.8357870x10-16 

-1. 5243436x10-17 

9 

1.4685126x10-1® 

6.7308017x10-30 

10 

1.7197282x10-21 

6.7792208x10-23 

11 

-4.3214861x10-34 

-1.0959706x10-33 

12 

-3. 5201050x10-37 

-7.8767464x10-30 

13 

0.0 

0.0 
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OHCsro'^invO  I  1  Of-lcNco^LOvD 


-7.5000000x10”^ 
3.7597711x10"^ 
-9.7111261x10"' 
4. 7984694x10"  j-J 
-7.4311015x10"^° 
4.2725291x10"^* 

-7.3171341x10"27 


-1.2500000x10"^ 

2.6855499x10"^ 

-4.4141460x10"° 

1.5995248x10"** 

-1.9408545x10"*^ 

8.6315720x10"*^ 

-1. 002898 5x10"2S 


-5.2083315x10"2 

4.1561584x10"^ 

-4.2136356x10"^ 

1.1054316x10"*^ 

-1.0269010x10"*° 

3.0250342x10"^^ 

0.0 


-3.9062497x10"^ 

-3.8625199x10"^ 

-3.8571553x10"^ 

1.0354265x10"?-^ 

-1.0007218x10"*° 

3.6710074x10"*° 

-3.5817803x10"°® 


-9.1145820x10"^ 

4.3293579x10"° 

-3.1061901x10"*® 

6.2916872x10"*^ 

-4.6890740x10"*® 

1.1506409x10"2° 

0.0 
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114.  f><»R*ll-H*l 
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IU6.  44  S40 

107.  V4P*14 

1L8.  DPS.f2*P*lJ/t2*P*21 
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Ot'5  ^lJ21  CaTANACH  solution  FUK  P  SUB  PHIME  AND  F  SUH 
AUG  if  9  67  14  49.4 
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FIGURE  3 


Stress  Intensity  Ratio  in  Tension  vs. A 
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FIGURE  4 


Stress  Intensity  Ratio  in  Bending  vs. 
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FIGURE  6 


Schematic  Diagram  of  Pressurization  System 


Electric  - 

Hydraulic  - 
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FIGURE  8 

Mesisured  Hoop  Strain  vs.  pressure 
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tance  From  center  of  crack 


FIGURE  12 
da/dn  vs .  K 
Shell  #10 


FIGURE  13 


da/dn  vs. 
Shell  #11 
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FIGURE  18 
da/dn  vs .  K 
Shell  #18 


(#/in 


FIGURE  19 


da/dn  vs. 


combined  Data 


Kr  (#/in.  3/2 


E 


FIGURE  20 


da/dn  vs .  k^. 
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fi72^ 


Measured  and  Theoretical  Strain  Ahead  of  Crack 


(  •ut/*ut^  ^  3 
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Na/2r 


Measured  and  Theoretical  strain  Ahead  of  crack 


Va/(2r) 


Measured  and  Theoretical  Strain  Ahead  of  crack 


p1a/2r 
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•  P  ia./2r 


f  Crack 

0 

a 

0 

*H 

ro 

P 

<u 

3 

P 

< 

•H 

P 

c 

P 

•H 

0 

nJ 

•H 

P 

CM 

•P 

w 

0 

Crl 

0 

cc: 

rM 

0 

ro 

U 

o 

0 

P 

H 

■H 

w 

Pm 

P 

o 

a 

u 

■H 

0 

P 

0 

0 

rC 

rd 

EH 

TJ 

iH 

w 

a 

p 

rd 

u 

(d 

0 

u 

0 

0 

0 

w 

{  ‘UT/'UT  W  ) 


100 


•  p  Va/(2r) 
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TABLE  I 


Ratio  of  Bending  and  Extension  Stress  intensity  Factors 


Lambda 

Number 
of  Terms 

Pi 

Pe 

0.2 

1 

-.0051019192 

1.0098910 

2 

-.0041031887 

1.0096026 

3 

-.0041033655 

1.0096165 

4 

-.0041033621 

1.0096165 

5 

-.0041033620 

1.0096165 

0.4 

1 

-.015236659 

1.0400153 

2 

-.011216673 

1.0365794 

3 

-.011225325 

1.0368121 

4 

-.011225110 

1.0368120 

5 

-.011225111 

1.0368120 

0.6 

1 

-.028057980 

1.0914843 

2 

-.018787392 

1.0770401 

3 

-.018868454 

1.0783070 

4 

-.018866012 

1.0783040 

5 

-.018866060 

1.0783044 

0.8 

1 

-.042834541 

1.1662898 

2 

-.025568961 

1.1262195 

3 

-.025954295 

1.1305623 

4 

-.025940092 

1.1305323 

5 

-.025940806 

1.1305365 

1.0 

1 

-.059402877 

1.2676320 

2 

-.030603158 

1.1788943 

3 

-.031862534 

1.1903994 

4 

-.031803183 

1.1902133 

5 

-.031808820 

1.1902435 

1.2 

1 

-.077912602 

1.4002587 

2 

-.033109799 

1.2299395 

3 

-.036358222 

1.2557787 

4 

-.036151880 

1.2549311 

5 

-.036181704 

1.2550822 

1.4 

1 

-.098740334 

1.5709013 

2 

-.032583775 

1.2749862 

3 

-.039741104 

1.3269937 

4 

-.039103933 

1.3238708 

5 

-.039224727 

1.3244645 
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Lambda 

Number 
of  Terms 

Pb 

Pe 

1.6 

1 

-.12246910 

1.7888593 

2 

-.028898230 

1.3108542 

3 

-.043134074 

1.4087889 

4 

-.041332584 

1.3988892 

5 

-.041742601 

1.4008612 

6 

-.041794437 

1.4009645 

7 

-.041820487 

1.4010472 

8 

-.041824244 

1.4010507 

• 

00 

1 

-.14989761 

2.0667491 

2 

-.022320944 

1.3355743 

3 

-.049097614 

1.5152988 

4 

-.044323031 

1.4869035 

5 

-.045575215 

1.4927968 

6 

-.045875596 

1.4936068 

7 

-.046006773 

1.4940039 

8 

-.046040097 

1.4940656 

o 

• 

CM 

1 

-.18205641 

2.4213584 

2 

-.013448964 

1.3482143 

3 

-.063462692 

1.6848021 

4 

-.051162538 

1.6069464 

5 

-.054819958 

1.6237831 

6 

-.056320724 

1.6288229 

7 

-.056919967 

1.6305416 

8 

-.057171478 

1.6312284 

2 

3 

4 

5 

6 

7 

8 


-.22020109 

-.003138382 

-.10295845 

-.069631625 

-.080504750 

-.087935212 

-.090599756 

-.092476415 


2.8743619 

1.3491155 

2.0391259 

1.8161060 

1.8654519 

1.8963699 

1.9035743 

1.9106433 


105 


TABLE  II 


MECHANICAL  PROPERTIES  OF  6063-T6  ALUMINUM 


Property 

Typical 

value 

Average  value 

Meridional 

for  Shell  #10 

Circumferential 

Strength  (psi) 

Tensile 

35,000 

29,260 

31,450 

Yield 

31,000 

27,790 

27,850 

Elongation 

in  2"  (%) 

12  -  18 

5.5 

E  (psi) 

o 

1— 1 

1.030x10^ 

.907x10'^ 

V 

1/3 

— 

TABLE  III 


DIMENSIONS  OF  6063-T6  ALUMINUM  SHELLS 


Nominal  Average  values 

Size  (in.) 

(in.)  Shell  shell  Shell  shell  Shell 
_ 9 _ 10 _ 12 _ n _ 15 

Diameter 

outside  8.00 


inside 

3.7880 

— 

3.7881 

3.7881 

3.7879 

variation 

-^.0008 

-.0006 

— 

+.0003 

-.0005 

+.0005 

-.0008 

+.0006 

-.0011 

wall  Thickness  .072 

— 

.0745 

— 

— 

— 

variation  — 

— 

-^.0006 

-.0010 

— 

— 

— 

106 


plat  to  shear  Fracture 
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TABLE  V 


crack  Growth  constants  of 
6063-T6  Aluminum 


Shell  0  A*  2(a3^+o(2)  rj^  log  Syx 

9  1.087  2.69xlO~J^  3.82  0.994  0.153 

10  1.066  6.99x10";^  3.99  0.980  0.284 

11  1.065  2.82x10”^®  3.46  0.970  0.303 

12  1.065  1.97x10"^®  3.54  0.965  0.339 

13  3.259  2.54x10"?-;  3.29  0.979  0.171 

14  1.077  1.16x10"?-;  3.37  0.994  0.127 

15  1.644  5.05x10"?-?  3.19  0.998  0.069 

18  1.055  4.21x10"^^  3.46  0.965  0.276 


_ig  .753  3.85 

da/dn  -  1.068x10  -^^(l+Q) 

A*  =  B(l+  0 


Kre  +  Krb 


lell 

u 

A' 

2{«i^-«2) 

^1 

log  Syx 

9 

1.087 

4.05x10"?-^ 

3.76 

0.994 

0.156 

10 

1.066 

1.00x10"?-^ 

3.93 

0.979 

0.284 

11 

1.065 

4.07x10"?-° 

3.40 

0.970 

0.304 

12 

1.065 

2.87x10"?-° 

3.48 

0.964 

0.341 

13 

3.259 

3.20x10"?; 

3.25 

0.978 

0.171 

14 

1.077 

1.61x10"?-; 

3.32 

0.994 

0.130 

15 

1.644 

6.86x10"?-? 

3.15 

0.998 

0.069 

18 

1.055 

6.09x10"^^ 

3.41 

0.965 

0.280 

da/dn  -  1.39x10“ 

^^(1-f  0)* 

729  3.80 
^r 

A‘  -  B(l+  0 ) 
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